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Abstract. We establish parabolicity and quadratic area growth for minimal 
surfaces-with-boundary contained in regions of R 3 which are within a sub- 
logarithmic factor of the exterior of a cone. Unlike previous work showing 
that these two properties hold for minimal surfaces-with-boundary contained 
between two catenoids, we do not make use of universal superharmonic func- 
tions. Instead, we use stochastic methods, which have the additional feature of 
giving a type of parabolicity in a more general context than Brownian motion 
on a minimal surface. 



1. Introduction 

We are interested in controlling the geometry and conformal structure of minimal 
surfaces-with-boundary in R 3 , under the assumption that they are contained in 
certain rotationally symmetric regions. In particular, let (x\, X2, X3) be standard 
Euclidean coordinates on R 3 and let r = \J x\ + x\. We consider minimal surfaces- 
with-boundary contained in sets of the form {|i3| < f(r) and r > C} for some 
positive, increasing, continuous function / and some positive constant C . We are 
interested in two possible properties of such surfaces, parabolicity and quadratic 
area growth. A surface-with-boundary is parabolic if the boundary is non-empty 
and any bounded harmonic function on the surface is determined by its boundary 
values. (This is not the standard notion of parabolicity for a Riemann surface, but 
it will be a convenient usage for us.) Equivalently, a surface is parabolic if Brownian 
motion on the surface, started from any interior point, hits the boundary (in finite 
time) almost surely. Now for positive p, let B p = {x\ + x\ + x\ < p 2 } be the ball 
of radius p. We say that a surface M has quadratic area growth if M n B p has 
area less than or equal to Cp 1 for sufficiently large p, for some positive constant C. 
(Note that this is an extrinsic notion of quadratic area growth.) 

We prove that both of these properties hold for minimal surfaces-with-boundary 
contained in regions that are within sub-logarithmic factors of the exterior of a cone 
(that is, f(r) as above is within a sub-logarithmic factor of being linear), under mild 
additional assumptions like stochastic completeness or properness of the immersion. 
These results are given in Theorems [3] and [H We prove both results using a sto- 
chastic approach, in contrast to the analytic approach via universal superharmonic 
functions (as described below). One additional feature of this approach is that 
parabolicity can be proven for a broader class of processes than Brownian motion 
on minimal surfaces, in a sense described in Theorem [2] 
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In the remainder of this section, we give additional background on this problem 
and make a few preliminary observations. We say that a surface M is contained 
between two catenoids if M is contained in the region {\xs\ < clogr and r > 2} for 
some c > 0. In [2], parabolicity and quadratic area growth were shown for prop- 
erly immersed minimal surfaces with compact boundary contained between two 
catenoids (see also Section 6 of the comprehensive survey [5]). This was done using 
the technique of universal superharmonic functions, where a universal superhar- 
monic function is a function on a subset of M 3 with the property that its restriction 
to any minimal surface is superharmonic. More specifically, the method of [2] is 
to find a universal superharmonic function on the region between two catenoids 
which is also proper and bounded from below. The existence of such a function 
immediately shows that any minimal surface of this type is parabolic. Further, 
using that the gradient and Laplacian of their universal superharmonic function 
are natural geometric quantities, they were able to deduce quadratic area growth. 
Bill Meeks has asked whether universal superharmonic functions which are proper 
and bounded from below can be found on the exterior of sufficiently large cones, 
mainly for the purpose of extending the parabolicity and quadratic area growth re- 
sults just mentioned. The results of the present paper bring us much closer (within 
sub-logarithmic factors) to the region conjectured by Meeks (of course, they leave 
open the question of the existence of universal superharmonic functions themselves 
on such regions). 

We also mention the complimentary project of finding regions (preferably as 
small as possible) where one can minimally (conformally) immerse a surface of 
arbitrary conformal structure. Recent progress in this direction can be found in pQ. 

The motivation for studying minimal surfaces-with-boundary restricted to this 
type of region comes in part from the fact that such surfaces arise as (representatives 
of) ends of complete minimal surfaces. In particular, the results of [2] on minimal 
surfaces-with-boundary contained between two catenoids, mentioned above, were 
used in the same paper in order to advance the classification of properly embedded 
minimal surfaces by controlling the behavior of certain ends. In light of this con- 
nection, as well as the fact that parabolicity and quadratic area growth are global 
properties, we are primarily interested in the large r behavior of the surfaces, and 
thus our results are formulated without much concern for the small r behavior of 
the surfaces. In this context, we note that parabolicity and quadratic area growth 
are preserved under homothetic rescaling (which also preserves minimality of the 
surface) and under removing compact subsets of the surface. This allows the re- 
sults proven below to be applied somewhat more broadly than the hypotheses might 
otherwise indicate. 

We make one more preliminary observation. Suppose that M is contained in the 
exterior of a cone, that is, M is contained in the region {(xi, x%, X3) '■ \xa\ < r} for 
some positive constant c. Let T p — {(x%, X2, £3) : r < p} for positive p (so that T p 
is a cylindrical region) . Then we see that M has quadratic area growth if and only 
if the area of M n T p is bounded by Cp 2 for large p, for some positive constant C. 
That is to say, quadratic area growth with respect to the distance from the origin 
is equivalent to quadratic area growth with respect to r for such surfaces. All of 
the surfaces we consider in this paper are contained in the exterior of a cone and 
thus satisfy this equivalence. In what follows, we will frequently use the condition 
on M D T p as our criterion for quadratic area growth without further comment. 
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2. Background on Markov chains 

A central tool for us will be the comparison of Brownian motion on M (or a more 
general martingale) to a Markov chain. In preparation for this, we briefly review 
some relevant facts about Markov chains, formulated for our case of interest. 

We take a subset of the non-negative integers, VLl = {L, L + 1,L + 2, . . .} for 
some non-negative integer L, as our state space. We are interested in discrete time 
Markov chains (Y n , n e {0, 1, . . .}) that evolve in the following way. If Y n = m > L, 
then Y n +i is m + 1 with probability p m and m — 1 with probability q m = 1 — p m . 
We stop the process upon hitting the boundary at L. So Y n is a time-homogenous 
nearest-neighbor random walk on f2^, stopped at L. We will always assume that 
p m > and q m > for all to > L, so that our chain is irreducible. Obviously, one 
could reduce this to the usual case of L = by a simple translation. However, this 
more general set-up will make the notation somewhat less burdensome later on. 

For to > L + 1 , we let 

, _ , QL + l gL+lg£+2 9L+1 •••qm 

Am — I H I r • • • H • 



PL+l PL+lPL+2 PL+l-'-Pm 

A well-known computation (see, for example, Theorems 5.3.7 and 5.3.8 of [5]) shows 
that A m is the reciprocal of the probability that the process, starting from L+ 1, 
hits to + 1 before L (we note that, under our assumptions, the process almost surely 
hits either L or to + 1 in finite time). Thus, the chain is parabolic, meaning that it 
almost surely hits L, if and only if A m — > oo as to — > oo. 

To give an example which will be relevant later, let q m = m/(2m + 1) and 
Pm = (m + 1)/(2to+ 1). Then 



q L +i---q m 1 , A 

= — and Am 

PL+\---Pm m hL 



Because the harmonic series is divergent, it follows that this chain is parabolic, and 
more generally, so is any chain with p m < {m + I) /(2m + 1) for large to. Further, 
we think of this chain as a borderline case, in the sense that the harmonic series is 
a borderline case for the divergence of series. 

In addition, we will need a result on an expected number of upcrossings. We 
assume that the chain is parabolic. Then if we start the process at Y = to, we let 
o\ be the first hitting time of to + 1, which may be infinite if the process reaches 
L first. Then we let u\ be the first hitting time of to after a\, 02 the first hitting 
time of to + 1 after a±, and so on. Note that if er^ is finite, then <7, is also finite, 
almost surely, because the chain is parabolic. If we let u m = sup{m : a m < 00} 
(with the convention that u m — if o\ =00), then u m is the number of upcrossings 
from to to to + 1. Because the chain is parabolic, u m is almost surely finite. A 
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straightforward adaptation of the computation that leads to A m shows that 
P (Y n hits L before m + 1\Y = m) 

^ I _j_ Pra PmPm-1 Pm ' ' ' PL+1 



q m q m q m -i q m ---qh+i. 

Further, by the Markov property, the number of upcrossings from m to m + 1 is a 
geometric random variable with the above as its parameter, and thus its expectation 
is 

TP, i v x 1 . Pm PmPm-1 . . Pm ■ ■ ■ PL+1 

E (u m \Y = m) = H 1 1 1 . 

q-m qmQm-1 Qm " " " 

3. Control of the radial process 

Having recalled some basic facts about Markov chains, we now discuss how they 
arise in the present work. We begin by introducing a broader class of processes, 
which we will use in our discussion of parabolicity. Consider an M 3 -martingale Z t 
which solves an (Ito) SDE of the form 

dZ t = a t dWf , where Wf is an R -Brownian motion, 

and where a t is an adapted process taking values in the set of 3 x 3 matrices 
conjugate under rotation to 

"1 0" 
10. 


That is, a t takes values in the set of projections onto planes. For simplicity, we will 
also assume that at is a continuous process. Intuitively, such a Zt is infinitesimally 
Brownian motion on the plane corresponding to a, where a is allowed to vary. 
One possibility is for a t to be projection onto the tangent plane of a minimal 
surface. More specifically, let M be a stochastically complete minimal surface-with- 
boundary. Stochastic completeness means that Brownian motion on M exists until 
it hits the boundary (that is, it doesn't explode), and we note that all properly 
immersed minimal surfaces-with-boundary are stochastically complete. Then Z t 
is Brownian motion on M, immersed into R 3 . (Here we are using the fact that 
such a Z t stays on the minimal surface, so a t is well-defined as projection onto 
the tangent plane at Z tl where self-intersections are handled in the usual way.) 
This is the motivating case. However, the point is that nothing in our approach to 
parabolicity will require the plane field corresponding to at to be integrable, so that 
Zt need not be restricted to a surface. In particular, in what follows it is appealing 
to consider the case when we choose cr t as a control in order to try to force Z t to 
remain in some subset of M? . As usual, we also wish to allow the possibility that 
our process Z t is stopped in finite time; in particular, Z t runs until some stopping 
time C, which may be infinite. In the minimal surface-with-boundary case, we will 
always take ( to be the hitting time of the boundary. 

For convenience, we call such a Z t a rank 2 martingale. The terminology is mo- 
tivated by the fact that at has rank 2, which means that the process infinitesimally 
evolves along a plane. Of course, we're assuming more than just that a t has rank 
2, corresponding to assuming that Xt is not merely infinitesimally evolving along a 
plane but is doing so "like Brownian motion." Nonetheless, making our terminology 
more specific at the expense of making it more cumbersome doesn't seem worth it. 
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Recall that we let r — \/x\ + x\ be the radial coordinate, in the sense of cylin- 
drical coordinates. We will assume that Z t is contained in the region {r > e L } 
for some non- negative integer L, and that it is stopped if it ever hits {r — e L }. 
For convenience, we will also assume that r(Z ) = e L+1 . We now discretize Z t (in 
particular, its radial component r t = r(Z t )) in the following way. Set = 0, and 
let (Ti be the first time logr t hits either L (in which case the process is stopped) 
or L + 2. Then we inductively define a n+ i for n > 1 to be the first time logr t hits 
either logr CTri + 1 or logr^ — 1. As usual, we allow a n to be oo if the event never 
occurs. 

Let rit be a unit vector in the kernel of at, and let 77,3^ be its a^-component. 
The notation is chosen because when Z t is Brownian motion on a minimal surface, 
n t is a unit normal to M (at the current point). Then a simple application of 
Ito's rule shows that a t = 1 — ttJ t is the infinitesimal drift of r\, meaning that 

r(Z t ) 2 — Jq at dt is a martingale. In the minimal surface case, this is just the fact 
that (Am/2)t 2 = 1 — n 2 , where 773 is the ^-component of the unit normal to M. 
Since 1 < a t < 2, it follows that the expectation of r\ grows at a rate between 
1 and 2. In particular, we see that logr t almost surely leaves any interval of the 
form (to + l,m — 1) for to > L + 1, unless the process is stopped first. Thus, 
almost surely, a n will be 00 if and only £ is finite and occurs before a n . So we have 
that X n — logr CTn is a nearest-neighbor random walk on {L, L + 1, . . .}, started at 
Xq = L+ 1 and stopped at the random time C (where this has the obvious meaning 
in terms of the underlying process Z t ). If the process is stopped before n, in the 
sense that £ < a n , we will say that X n does not exist. Note that X n is not, in 
general, Markov. This is true even in the case when Z t is Brownian motion on a 
minimal surface-with-boundary M, since the transition probabilities (as well as £) 
depend on where Z an is (on M). 

We now further assume that Z t is contained in the region 

{M</(r)} 

for some continuous, non-negative, monotone non-decreasing function / : [e L , 00) 
[0,oo). Our goal is to find a Markov chain with which to compare X n . The first 
step is to find upper bounds on V (X n+ i =m + l\X n = to) that depend only on 
{X n = to} and not on anything else in (here T% is the c-algebra generated 
by Z t up to time a n ). 

Recall the following basic estimate for minimal surfaces (see Lemma 2.3 of [2]) 

Lemma 1. Let M be a minimal surface-with-boundary and assume that r 7^ on 
M. Then |A M k>gr| < |VMX 3 | 2 /r 2 (on the interior of M). 

The key observation for us is that this generalizes immediately to rank 2 mar- 
tingales. We let Pt and j t be the infinitesimal drifts of logrt and x\ t — x^(Z t ), 
respectively, analogously to our earlier introduction of a t ■ Then the generalization 
is 

(1) |/3 4 | < — J assuming that r t 7^ 0. 

For completeness, we now take a moment to derive this inequality. Because at is a 
projection, of = at- Thus, in terms of a unit vector (77,1,77,2,713) in the kernel of a t , 
the diffusion matrix associated to Z t is (and where we feel free to drop the t's in 
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the subscripts to simplify the notation) 



1 - n\ 
—n\n<2, 
-nxn 3 



—n\ni —n\ri3 
1 - n| -n 2 n 3 



-n 2 n 3 1 

The Hessian of log r (where r ^ 0) is given by 



d 2 



l 



dx\ ^°& r r 2 



2x{ 

t-4 



if 



dx\dx2 
dx\ 



log r — 



d 2 



dxidx\ 



logr = 



2x\X2 



Then Ito's rule gives 



1 

2^ 



logr = 



2xj_ 



1 



2x 2 



and for everything else. 



2x 2 



Similarly, we compute that jt = 1 — Tt§ = n\ + n|. Because the desired inequality 
is invariant under rotation around the a^-axis, it is enough to prove it under the 
assumption that x\ = r and x 2 = at the current point. Making this assumption, 
we get the simplification fit = (nl—n^/Zr 2 . Comparing this with 7t/2r 2 establishes 
the desired inequality. In the case when Z t is Brownian motion on M, we have that 
fit = (Am/2) logr and "f t = |V ' mXz\ 2 , where the right-hand sides are evaluated 
at Zt, of course. Thus the situation reduces to that of Lemma HJ interpreted 
probabilistically. 

Next, our assumptions on / (and the definition of the a n ) imply that 



E 



7t dt 
E [x 2 (Z a 



{X n = m} n T'i 



<(/(e^)) 



+lAC. 
2 



~x 2 (Z an )\{X n = m}n^J 



Note that this bound is independent of everything in J 7 ^ except {X n = to}. We 
use this estimate along with Equation (Q} and Ito's rule to see that 



E [log 



\{x n = m}n^ n ] < 



(f(e m+1 )Y 



2e 



2m— 2 



which again is independent of everything in T% except {X n = m}. Recall that, 
given X n = m, X n +\ = logr CTrl+1 is either m+1 or m — 1, assuming er„+i < oo (that 
is, the process isn't stopped first), and if <7 n +i = o° then logr^ G [m — l,m + 1]. 
Thus, we see that 



i\x n 



, 1 



(/ (e m+1 )r 

4 e 2m-2 



independent of anything else in J 7 ^ . In particular, we have found the desired 
uniform estimate for P (X n +i — m + l\X n = to). 

We wish to use this bound as the transition probability p m for our comparison 
chain. This means that it must be less than one (recall that we want our chain to 
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P'»=o+ , o2m - 2 for all to > L + 1. 



be irreducible). Thus, we will now also assume that L and f(r) are such that 

(f(e m+1 )) 2 1 
(2) l2 m -2 <2 ^ a \\m>L + l. 

(Note that this is done mainly for convenience.) We set 

1 , (f(e m+1 )f 

2 4e 2 

Next, we assume that the probability space on which Zt is defined is rich enough 
to support a countable sequence U\, U-x, ■ ■ ■ of i.i.d. random variables, distributed 
uniformly on the interval [0,1], and all independent of Zt (this can always be 
accomplished; for example, by taking the product of the original probability space 
with the unit interval with Lebesgue measure). We now determine the process Y n 
as follows. Let Y = L + 1, and define Y n for n > 1 inductively as follows. If 
Y n _i = m > X n -\ or if X n _i doesn't exist (equivalently, Q < cr n _i), then 



Yn = 

If Y n -i = m = X n -i, then let 

ifi m = P \X n 



?7i + l if U n € [0,p m ], 
771-1 if U n e (pm, !]• 



Note that ip m is bounded from above by p m (so ip m < p m ), and let 

{ 777 + 1 if X n = 771 + 1 , 
m + 1 if X n ^ to + 1 and U„ < (p m - ip m ) / (1 - <p m ), 
771-1 if X n ^ to + 1 and U n > (p m - ip m ) / (1 - ^ m ). 

Here X n ^ m + 1 should be thought of as short-hand for X„ = m — 1 or X n doesn't 
exist. Finally, we stop Y n if it hits L. Recall that we're also assuming that Z t is 
stopped if it hits {r = e L }. 

A simple computation shows that P [Y n = m + l|y n _i = m] = p m , independent 
of anything else in _i an( i °f • • • > U n —\. It follows that y„ is a Markov chain 
with transition probabilities p m . In addition, by construction Y n dominates X n in 
the sense that Y n > X n for all n such that a n < C, almost surely. So Y^, is the 
desired comparison chain. 



4. PARABOLICITY 

Because the conditions for minimal surfaces-with-boundary we are ultimately 
interested in, namely parabolicity and quadratic area growth, aren't affected by 
homothetically rescaling (as mentioned in the introduction), we will assume that 
M is contained in the region {r > e L } for some non-negative integer L (instead 
of the more general {r > c} for positive c), and that any intersection of M and 
{r = e L } is contained in the boundary of M. 

Let 

CT 

fi(r) = —j^^^=^^^= for some c > 0, 
v /logrlog(logr) 

and let L be chosen so that Equation ([2]) is satisfied with f — fi- Note that a 
simple computation shows that such an L always exists. (We should perhaps use 
Li here, analogously to /i, but that would make the subscripts too unwieldy.) 
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The point of working with rank 2 martingales in the last section is that we have 
the following more general version of parabolicity. 

Theorem 2. Let fi(r) and L be as in the preceding paragraph. Suppose that Z t is 
a rank 2 martingale as described in Section^ (in particular, we continue to assume, 
for convenience, that r(Zo) = e L ). If Zt is contained in the region 

{r > e L and \x 3 \ < fi{r)} C R 3 , 

then Zt almost surely has finite lifetime. 

Proof: We recall X n , the discretization of \ogr(Z t ) introduced above. Then with 
/ = fi, our earlier computations show that we have a comparison Markov chain 
Y n , started at Yo = L + 1, with (after a bit of simplification) 

1 c¥ 1 
Pm ~ 2 + ~ ' (m + 1) log(m + 1) ' 
Then a little more algebra shows that 

Prn < if and only if -5— r (m + 1) log(m + 1) > 4m + 2. 

2m + 1 ere 4 

This last inequality holds for all sufficiently large m. Thus p m < (m+ l)/(2m+ f) 
for large m, and it follows from the discussion in Section [5] that Y n is parabolic, 
that is, Y n almost surely hits L. 

Recall that Y n dominates X n until £. Then because Zt is contained in the region 
{r > e L } and we stop Zt when log r(Zt) = L, if not before, Zt is stopped no later 
than the time when Y n hits L, where this has the obvious meaning in terms of the 
stopping times o~ n . In particular, £ is almost surely finite. □ 

Returning to minimal surfaces, recall that asserting that a surface- with-boundary 
is parabolic implies that the boundary is nonempty. Again, we are interested in the 
asymptotic behavior of minimal surfaces, so the following theorem is formulated 
without much concern for the behavior of the minimal surface for "small r." 

Theorem 3. Let /i(r) and L be as above. Suppose that M is a stochastically 
complete minimal surface-with-boundary contained in the region 

{r>e L and \x 3 \ < /i(r)} C R 3 , 

with the assumption that that any intersection of M and {r — e L } is contained in 
the boundary of M . Then M is parabolic. 

Proof: We consider Brownian motion B t on M started at some interior point 
with r = e L+1 , which we can assume exists after homothetically rescaling (though 
we might have to change c and L after rescaling) and using that M is stochastically 
complete. We let £ be the first hitting time of the boundary and stop B t at £. Recall 
that, in order to show that M is parabolic, it is enough to show that Brownian 
motion started at any interior point hits the boundary in finite time, that is, Q is 
almost surely finite. Now B t = Z t satisfies the hypotheses of Theorem [21 and we 
conclude that B t almost surely hits the boundary in finite time. □ 

The same arguments, along with more involved computations, could presum- 
ably be pushed to give a slightly weaker hypotheses, meaning replacing f\{r) with 
something larger by a sub-logarithmic factor. However, these arguments, in their 
present form, break down if we try to extend them to the case of surfaces (or rank 
2 martingales) contained in the exterior of a cone (that is, the case when fi(r) 
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is replaced by cr for some c > 0). The reason is that such a region is preserved 
under homothetic rescaling, and thus the p m (as defined above) don't decay to 1/2 
as m — > oo. In another direction, one could estimate quantities associated with 
the parabolicity of M, for example, the distribution of the maximum of r along 
Brownian paths. However, we have not pursued these directions here because we 
are unaware of any applications of this kind of result. 

The reason for framing the discussion in terms of rank 2 martingales is that, de- 
pending on your point of view, it provides an appealing "explanation" for parabol- 
icity, and clarifies the hypotheses. In particular, it is interesting that parabolicity 
doesn't rely on having a surface at all, or equivalently, on the plane field determining 
the evolution of the rank 2 martingale being integrable. 

5. Area growth 

While it is common to study parabolicity using Brownian motion, it is perhaps 
more interesting that it also makes a fairly natural tool, at least in this context, to 
study area growth. The generalization of quadratic area growth to rank 2 martin- 
gales would be a quadratic estimate on occupation times, as the proof the following 
theorem makes clear. However, working with rank 2 martingales doesn't seem 
worth the added complexity in the present context. Thus, for the remainder of the 
section we restrict our attention to the case when our rank 2 martingale is Brown- 
ian motion on a stochastically complete minimal surface- with-boundary, and adopt 
notation accordingly. 

Let 

CT 

f 2 {r) = —==- — for some c> 0, 

Vlog r log (log r) 

and let L be chosen so that Equation @ is satisfied with / = f%. Again, a simple 
computation shows that such an L always exists. Note that for a minimal surface 
satisfying the assumptions of the following theorem, quadratic area growth is equiv- 
alent to quadratic area growth with respect to r, as discussed in the introduction. 

Theorem 4. Let f2(r) and L be as in the preceding paragraph. Suppose that M is 
a properly immersed minimal surface with compact, non-empty boundary, contained 
in the region 

{r > e L and \x 3 \ < / 2 (r)} C E 3 

with the assumption that any intersection of M and {r = e L } is contained in the 
boundary of M. Then M has quadratic area growth. 

Proof: Removing a compact subset of M doesn't affect the result, nor does 
homothetic rescaling (though we might have to change c and L after rescaling). 
Thus, without loss of generality, we can assume that the boundary of M is equal to 
M n {r = e L } and that M n {r = e L+1 } is a smooth, compact curve in M, which 
we denote 7. Then 7 divides M into two connected components (which intersect 
only at 7), the unbounded piece M u = M fl {r > e L+1 } with boundary 7, and the 
compact piece M c = M n {e L < r < e L+1 } with two boundary components, 7 and 
dM (the original boundary of M). Also, since dM = M n {r = e L }, we see that 
Brownian motion is only stopped when it hits the level {r = e L }. 

Note that for large r, fa{r) < fi(r), regardless of what constants c are used in 
the two functions. So applying Theorem[3]and using that M is properly immersed, 
we see that M is parabolic. It follows that M u is also parabolic. 
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We consider the function h on M, determined as follows. We let h be harmonic 
on M c with boundary values 1 on 7 and on dM. Because M c is compact, these 
boundary values uniquely determine a function. On M u , we let h be the unique 
bounded harmonic function with boundary values 1 on 7 (which is well-defined 
because M u is parabolic); thus h is identically 1 on M u . Finally, we let h be 1 on 
7. Thus, h is a continuous, bounded function on M and is harmonic on M \ 7. 

Observe that ^Ah is measure- valued, supported on 7, and has finite mass; let 
that mass be I /a. Then |A(a/i) is a probability measure, which we denote p,, 
supported on 7. It follows that ah is the Green's function associated to p, and 
thus ah is also the density of the expected occupation time of Brownian motion on 
M started from p and stopped at dM. (See [4] for background on the relationship 
between Brownian motion and potential theory on manifolds.) Since we know that 
M c has finite area and ah is bounded there and that ah is identically equal to a 
on M u , it follows that the expected occupation time of M n {r < p} is asymptotic 
to a times the area of M D {r < p] (as p ~ > 00). This is the connection between 
Brownian motion and area growth that we need. In particular, suppose that, for 
Brownian motion started at some point on 7, we show that the expected occupation 
time of M n {r < p} is bounded from above by Cp 2 for some positive constant C. 
If this holds with C independent of which point of 7 our Brownian motion starts 
from, then the estimate will hold for Brownian motion started at 7 and we will 
have shown that M has quadratic area growth. 

To do this, let p = e k for some integer k > L + 1. We start be observing that 
Brownian motion started at a point with r = e L+1 (that is, a point on 7) almost 
surely hits either r = e L (in which case the process is stopped) or r — e k+1 (that 
is, ep) in finite time. Next, we need to estimate how long this takes. Let 9j be first 
hitting time of r — e 3 for any integer j; then 



E 



2 l„_ _ „L+l 

+ 1 



e 2 ( fc+1 )p (6 k+1 < 6 L ) + e 2L P (0 L < 6 k+1 ) 



(independent of where on {r = e L+1 } the Brownian motion starts). On the other 
hand, from 1 < (Ajn/2)r 2 < 2 and Ito's formula, we see that 



E 



re,AO h+1 \ro = e L+1 ] > e 2 ^ + E [6l A 9 k+1 ] 



(again independent of where the Brownian motion starts). It follows that 

E [6 L A 9 k+1 \r = e L+1 ] < e 2 ( fe+1 ) - e 2 ( L+1 >. 

This is the desired estimate on how long it takes for Brownian motion to leave 
r G (e L ,e k+1 ); note that it is quadratic in the upper endpoint e k+1 = ep and 
independent of where on 7 = {r = e L+1 } the Brownian motion starts. 

Continuing, we see that if 9l < Ok+i then the process is stopped at 8l and 
the only contribution to the occupation time of M n {r < p} is 9^. (Because 
we're stopping at r = e k+1 rather than e k , this is just an upper bound on the 
contribution; the reason for this set-up will be clear below.) If 8 k +i < #l then the 
process continues, and we must also estimate the contribution to the occupation 
time after 9 k +i- Intuitively, we do this as follows. Once the process hits r — e fe+1 , 
we stop counting it toward the occupation time, until the next time it hits r = e k 
(which almost surely happens). Then we start counting the time again, until the 
next time it hits either e L or e fe+1 , at which point either the process is stopped or we 
perform another round of the same procedure. Note that, almost surely, the process 
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eventually hits r = e L and is stopped. This procedure counts all of the occupation 
time of M n {r < p} and part of the occupation time of M D {p < r < e k+1 }, 
but it certainly provides an upper bound on the occupation time of M n {r < p}. 
Further, for the purpose of establishing quadratic area growth, this is no real loss, 
since e L+1 = e • p. 

We now formalize the above procedure. Let £o = #fc+i be the first hitting time of 
r = e k+1 , and let £o be the first hitting time of r = e k after £1. Then we inductively 
define to be the first hitting time of r = e k+1 after £ n _i and to be the first 
hitting time of r = e k after £„. By parabolicity, £ n is almost surely finite if £ n is. 
Let Uk — sup{n > : £ n < 00} be the number of upcrossings from r = e k to e fe+1 ; 
again by parabolicity it is almost surely finite. An argument completely analogous 
to the above shows that, starting from the expected amount of time spent 

until the process leaves r € (e L , e k+1 ) satisfies the estimate 

>l A £„) - |„_i| |„_i < 00] < e 2(fc+1) - e 2k . 

(We note that this estimate is uniform over where on {r = e fc } the process is at 
£n-i-) Thus, the contribution of each upcrossing to the expected occupation time 
is bounded by the expression on the right. 

Summarizing our progress so far, we have that the expected occupation time 
of M n {r < p}, where p = e k , for Brownian motion started at any point of 7 is 
bounded from above by 



eW J _ e ^ + l, + p {f)k+i < H) _ e ^ E [Ukl0k+i < 0l] 

Observe that e 2 ^ k+1 ^ — e 2 ( i + 1 ) and e 2 ( k+1 ) — e 2k are both quadratic in p — e k . 
Thus, in order to show that this expression grows quadratically (in p), it is enough 
to show that 

P(0 fe +i <9 L )E[U k \6 k+1 <9 L ] 

is bounded as k — > 00. Again, we recall X n , the discretization of log rt introduced 
above. We see that 9k+i corresponds to the first time X n hits k + 1 (where this 
correspondence is understood in the natural way in terms of the stopping times a n ) 
and Ojj corresponds to the first time X n hits L (and is stopped). Also, note that 
C = Ojj, so X n is a nearest-neighbor random walk on {L, L + 1, . . .} stopped when, 
and only when, it hits L. Thus we will have no need to refer explicitly to £, and 
the discussion in Section [3] simplifies accordingly in this case. (This simplification 
is made possible because here, unlike in Theorem [3J we are assuming that M is 
properly immersed with compact boundary, as discussed at the beginning of the 
proof.) 

We see that P(8k+i < 6l) can be re-expressed as the probability that X n hits 
k+ 1 before L. Further, letting f — f 2 , our earlier computations show that we have 
a comparison Markov chain Y n with 

_1 cV 1 

Pm ~ 2 ~T~ ' (m + 1) (log(m + l)) 2 ' 

Thus P (9k+i < 6l) is bounded from above by the probability that Y n , started from 
L+l, hits k+l before L, independent of where on 7 the underlying Brownian motion 
in M starts. Further, this probability was given in Section [5] as the reciprocal of 
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Ak, and recalling that gives 
P(6 k+1 <6 L ) < (l 



gL+l qL + lQL+2 gL+1 ' • • gfc 

PL+1 PL+lPL+2 PL+i---Pk 



where above and q m = 1 — p m , as usual. Also, observe that both X. n and 

F„ are parabolic, which follows from the proof of Theorem [31 

Next, we wish to estimate E \Uk\6k+i < 9l] m terms of Y n as well, and note that 
such an estimate in terms of Y n will hold for X n independent of where on 7 the un- 
derlying Brownian motion starts. The point is that, even though upcrossings of X n 
might occur at different times than upcrossings of Y n , parabolicity and the Markov 
property (of Y n ) imply that E [Uk\dk+i < 6l] is no larger than E (v,k\Y Q — k). This 
is more or less obvious, but we provide a bit more detail for the sake of clarity. The 
main idea is that we can feel free to pause either X n or Y n in order to allow the 
other process to "catch up" without changing the distribution of the total number 
of upcrossings of either process. First, because we're conditioning on 9k+i < &l, 
we know that X n will hit k + 1 and then by parabolicity will come back to k. We 
wait until the happens to start Y n (which we start from k). Then Y dominates 
X, so if X has an upcrossing (that is, if X hits k + 1 again before hitting L), so 
does Y. Supposing that X does have an upcrossing, it might come back down to 
k before Y does. However, once X comes back down to k, we can "pause" it there 
and wait for Y to also come back down to k, which it will by parabolicity. With 
both processes at k, we simply repeat the above procedure. This shows that every 
upcrossing of X corresponds to an upcrossing of Y, and thus the desired inequality 
for the expected numbers of upcrossings holds. 

Now using our computation of E (tife|io = k) from Section [21 we see that 



h+i <e L )E[u k \0 k+1 <e L ] < 



1 



1 Pk PkPk-i ^ Pk ■ ■ ■ Pl+i 

qk qkqk-i q k ---q L +i 

-I -1 

9l+i . 9L+19L+2 . . <Zl+i ' ' • qk 



PL + 1 PL + lPL+2 PL+l"'Pk 



Recall that, in order to show that the expected occupation time of M D {r < p}, 
where p = e k , grows quadratically in p, it is enough to show that the above quantity 
stays bounded as k — > 00. The right-hand side is a fraction (though not displayed 
as such), and because p m > q m the terms in the numerator are increasing (from left 
to right) while the terms in the denominator are decreasing. Since both numerator 
and denominator have the same number of terms (for a given k), we see that in 
order to show that this fraction is bounded as k — > 00, it is sufficient to show that 
the ratio of the right-most terms, that is, 

Pk---PL+i\ fqL+i---qk\ fPk---PL+i^ 2 



qk---qh+ij / \PL+i---Pk/ \q k ---qL+i 

is bounded as k — > 00. Taking the logarithm of the right-hand side, we see that it 
is enough to show that Yl m =L+i l°g(Pm/?m) is bounded as k — > 00 (recall that all 
terms in the sum are positive since p m > qm)- To do this, we will show that 

I'M log {^) < . 

1 q m J (m + 1) (log(m + l)) 2 
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for some c > 0, for large m (so that Yli=L+i ^°s(Pm/<lm) is bounded by the usual 
comparison test for sums). 

We have (using the above formula for p m ) that 

log ( — ] = log 




For large m, both p m and q m become arbitrarily close to 1/2. Thus we can use an 
estimate based on the first-order Taylor expansion of the logarithm, namely that 
| log(a; + 1/2) — log(l/2)| < 3|a;| for x sufficiently close to 0. This implies that, for 
large m, 



QmJ V 2 / 4 (m + 1) (log(m + 1)) 

, , n 3c 2 e 4 
log 



2 



,27 ' 4 (m + 1) (log(m + l)) 2 

_ 3cV 1 

2 (m + l)(log(m + l)) 2 ' 

This establishes (*). Combining all of the above, we have proven the desired es- 
timates on Brownian motion uniformly with respect to where on 7 the Brownian 
motion starts, and we can now conclude that there exists a positive constant C 
such that the area of M <~){r < p}, when p = e k , is bounded from above by Cp 2 for 
sufficiently large p. 

Finally, to extend this result to all p (and not just those that are an integer 
power of e, as we have been assuming to make use of our discretization X n of the 
radial process), suppose that p £ (e k ~ 1 , e k ]. Then since the area is monotone in p, 
we have 

Area(M n {r < p}) < Ce 2k = CeV^" 1 ) < (C*e 2 ) p 2 . 

So increasing C by a factor of e 2 gives the estimate for general p, completing the 
proof. □ 

Note that, having proven the above result, a standard argument can be used 
to extend it somewhat. In particular, the monotonicity formula for area (see [6]) 
shows that 

Area(M n {r < p}) 

2 

np z 

converges to a finite limit as p — > 00, which wc denote n(M). Then because M is 
contained in 

{r>e L and \x 3 \ < f 2 {r)} , 

geometric measure theory allows us to conclude that the homothetic shrinkings of 
M converge to a locally finite minimal varifold supported on the (xi, a;2)-plane. 
This limit is thus an integer multiple of the plane, and we conclude that n(M) is a 
positive integer. (See either the end of Section 2 in [2] or Section 6.4 of [5] for more 
details of the argument.) When M is an appropriate representative of an end of a 
minimal surface, n(M) is generally called the multiplicity of the end. 
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Recall that if the area of a surface grows quadratically with respect to the intrin- 
sic distance, then the surface is parabolic. In light of this, it is not surprising that 
these methods give parabolicity under a strictly weaker condition than is needed 
for quadratic area growth (with respect to r, or equivalently the exterior distance 
from the origin). Also, as was the case for Theorem [3l these computations could 
likely be extended to allow replacing fz(r) by something slightly more lenient, but 
they cannot be simply extended to allow /^(O to be replaced by /i(r), much less 
to allow M to be contained in the complement of a cone. Further, in various cases 
one could give explicit estimates for the constants appearing in the proof in order 
to derive an estimate for C, and thus also for n(M). In particular, a depends only 
on the behavior of a compact subset of M, such as M PI {r < e L+2 }. Moreover, 
P (9k+i < 9l) E [Uk\0k+i < 6l] can be estimated in terms of the p m , which depend 
only on fi- Replacing f% with a smaller / would decrease the p m and thus also the 
upper bound on n(M). However, as was the case previously, we have not pursued 
this because we are unaware of any interest in such an estimate. 
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